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Abstract: A survey of atomic binding energies used by general purpose Monte Carlo systems is 
reported. Various compilations of these parameters have been evaluated; their accuracy is estimated 
with respect to experimental data. Their effects on physics quantities relevant to Monte Carlo par- 
ticle transport are highlighted: X-ray fluorescence emission, electron and proton ionization cross 
sections, and Doppler broadening in Compton scattering. The effects due to different binding ener- 
gies are quantified with respect to experimental data. The results of the analysis provide quantitative 
ground for the selection of binding energies to optimize the accuracy of Monte Carlo simulation in 
experimental use cases. Recommendations on software design dealing with these parameters and 
on the improvement of data libraries for Monte Carlo simulation are discussed. 
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1. Introduction 

The simulation of particle interactions in matter involves a number of atomic physics parameters, 
whose values affect physics models applied to particle transport and experimental observables cal- 
culated by the simulation. Despite the fundamental character of these parameters, a consensus has 
not always been achieved about their values, and different Monte Carlo codes use different sets of 
parameters. 

Atomic parameters are especially relevant to simulation scenarios that are sensitive to detailed 
modeling of the properties of the interacting medium. Examples include the generation of char- 
acteristic lines resulting from X-ray fluorescence or Auger electron emission, and precision sim- 
ulation studies, such as microdosimetry, that involve the description of particle interactions with 
matter down to energies comparable with the scale of atomic binding energies. 
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Simulation in these domains has been for an extended time the object of specialized Monte 
Carlo codes; some general purpose Monte Carlo systems have devoted attention to these areas, in- 
troducing functionality for the simulation of fluorescence, PIXE (Particle Induced X-ray Emission) 
and microdosimetry. In this context, emphasis has been placed on the development and vaUdation 
of the physics models implemented in the simulation systems, while relatively limited effort has 
been invested into verifying the adequacy of the atomic parameters used by general purpose Monte 
Carlo codes with regard to the requirements of new application domains. 

This paper surveys atomic binding energies used by well known Monte Carlo systems, includ- 
ing EGS [1], EGSnrc [2], Geant4 [3,4], ITS (Integrated Tiger Series) [5], MCNP/MCNPX [6,7] 
and Penelope [8], and by some speciaUzed physics codes. These software systems use a variety of 
compilations of binding energies, which are derived from experimental data or theoretical calcula- 
tions; this paper investigates their accuracy and their effects on simulations. 

To date, to the best of our knowledge, no comprehensive review of the binding energies used 
by major Monte Carlo codes has been documented in the literature yet. The most recent review of 
atomic binding energies is included in the Handbook of X-Ray Data [9]: it is Umited to comparing 
two compilations of binding energies [10, 11], derived from experimental data, with respect to a 
common reference [12]; the comparison rests on the visual appraisal of plots. Out of the three 
compilations of binding energies examined, only two are used by major Monte Carlo systems. We 
are not aware either of previously published objective estimates, based on statistical methods, of 
the compatibility with experiment of the binding energies used in Monte Carlo simulation, and of 
simulation observables that depend on them. 

A small subset of preliminary results of this study were summarized in a conference paper [13]. 

2. Compilations of electron binding energies 

The binding energies considered in this study concern neutral atoms in their ground state; several 
compilations of their values, of experimental and theoretical origin, are available in the Uterature. 

Compilations based on experimental data are the result of the appUcation of selection, evalu- 
ation, manipulations (like interpolation and extrapolation) and semi-empirical criteria to available 
experimental measurements to produce a set of reference values covering the whole periodic system 
of the elements and the complete atomic structure of each element. 

Most of the collections of electron binding energies based on experimental data derive from 
a review pubUshed by Bearden and Burr in 1967 [10]. Later compilations introduced further re- 
finements in the evaluation of experimental data and the calculation of binding energies for which 
no measurements were available; they also accounted for new data taken after the pubUcation of 
Bearden and Burr's review. 

Experimental atomic binding energies can be affected by various sources of systematic effects; 
they originate not only from the use of different experimental techniques in the measurements, but 
also from physical effects: for instance, binding energies of elements in the soUd state are different 
from those of free atoms, and binding energy measurements can be affected by the chemical state 
of a solid. 

The first attempt to calculate electron binding energies was reported by Slater [14]; since 
then, various relativistic computations of neutral atom binding energies have been performed [15]. 
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They exploit methods based on a Dirac-Hartree-Slater model, with corrections for QED (quantum 
electrodynamics) effects and the nuclear charge distribution. 

2.1 Selected compilations 

This paper evaluates a selection of binding energy compilations, which are used by general purpose 
simulation systems and some representative speciaUzed codes: 

• the compilation by Bearden and Burr [10], 

• the compilation by Carlson [16], 

• the tabulation included in Evaluated Atomic Data Library (EADL) [17], 

• the compilation assembled by Sevier in 1979 [12], 

• the compilations included in the seventh and eighth editions of the Table of Isotopes [18, 19], 
respectively pubhshed in 1978 and 1996, 

• the compilation by Williams included in the X-ray Data Booklet [20] and in the CRC Hand- 
book of Chemistry and Physics [21]. 

Bearden and Burr performed a comprehensive evaluation of experimental X-ray wavelength 
data; the techniques they used to establish a consistent energy scale and to deal with elements with 
multiple or missing measurements are documented in [10]. This compilation has been the basis 
for several other ones pubhshed in the following years and is still used in some physics software 
systems. 

Carlson's compilation reproduces the one by Lotz [22] with a few modifications and exten- 
sions, that concern the data for krypton and xenon, the binding energies of elements with atomic 
number greater than 94 and the P shell data of elements with atomic number between 87 and 95. 
The compilation covers atomic numbers from 1 to 106; values are given for free atoms and are 
referenced to the vacuum potential. 

The compilation by Lotz is based on Bearden and Burr's evaluated data, complemented by 
other experimental measurements. The tabulated binding energies were determined according to 
empirical criteria, interpolation and extrapolation of available data. Since values are listed for free 
atoms, the work function was taken into account in converting experimental binding energies for 
sohds. According to [22], the uncertainties of the tabulated values are at most 2 eV for elements 
with atomic number up to 92; larger uncertainties, in some cases greater than 10 eV, are reported 
for heavier elements. 

The binding energies collected in the seventh edition of the Table of Isotopes (identified in the 
following as Tol 1978) were taken from Shirley et al. [23] for elements with atomic number up to 
30, and from a compilation of experimental data by the Uppsala Group [24] for heavier elements. 
The tabulated binding energies derive mainly from photoelectron spectroscopic measurements; 
data were taken from Bearden and Burr's compilation in cases where experimental photoelectron 
measurements were not available. Interpolation and extrapolation techniques were used to comple- 
ment experimental data. The data are listed with reference to the Fermi level and concern elements 
with atomic number from 1 to 104. Uncertainties are reported as about 0.1 eV for light elements 
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and 1-2 eV for most elements with atomic number greater than 30; uncertainties approaching 100 
eV are mentioned for transuranic elements. Shifts of the order of 10 eV in the binding energy of 
non- valence shells can result from changes in the chemical state of the medium [23]. 

The binding energies collected in the eighth edition of the Table of Isotopes (identified in the 
following as Tol 1996) were taken from the compilation by Larkins [25]. Binding energies are 
reported for soUd systems referenced to the Fermi level, except those for noble gases, CI and Br, 
which are for vapor phase systems referenced to the vacuum level. Uncertainties may be as large 
as 10-20 eV for the inner orbitals in the high-Z elements, and changes in chemical state can lead to 
substantial shifts in the binding energies of non- valence shells [23]. 

The binding energies tabulated by Larkins are based on Sevier's 1972 compilation [26] for 
elements with atomic number up to 83 and on the compilation by Porter and Freedman [27] for 
heavier elements; with respect to these references, Larkins includes some updated values for Ar, 
Ge, As, Se, Xe and Hg. Sevier's 1972 tabulations were mainly an update to Bearden and Burr's 
ones to include more recent measurements; a further extension was published by Sevier in 1979 
[12]. Porter and Freedman combined a theoretical approach and experimental measurements to 
interpolate data for heavy elements. 

The eighth edition of the Table of Isotopes also includes a Ust of ionization energies of the 
elements (concerning the least bound electron), which reflects the data available from NIST (United 
States National Institute of Standards and Technology) [28]; these values differ in some cases from 
those in the tabulation of electron binding energies in the same volume. 

Williams' compilation is based on Bearden and Burr's data; some values are taken from [33] 
with additional corrections, and some from [34]. The energies are given relative to the vacuum 
level for the rare gases and for H, N, O, F and CI, relative to the Fermi level for metals and relative 
to the top of the valence bands for semiconductors. The tabulations concern elements with atomic 
number between 1 and 92. 

The atomic subshell parameters collected in EADL are derived from theoretical calculations 
by Scofield [29,30]; besides these two references, EADL documentation cites a "private commu- 
nication" by Scofield, dated 1988, as a source of the data. Due to the scarcity of documentation 
about the origin of the binding energies listed in EADL, it is difficult to ascertain how they were 
calculated, and what assumptions and approximations may be underlying. Binding energy values, 
although not for all elements and shells of the periodic system, are reported in some publications 
by Scofield [31,32]; those in [32] appear consistent with EADL tabulations. EADL data concern 
isolated, neutral atoms with atomic number up to 100. 

2.2 Binding energies used by physics software systems 

General purpose Monte Carlo systems and specialized codes use a variety of binding energy com- 
pilations. 

EGS5 uses the binding energies tabulated in the 1996 edition of the Table of Isotopes, while 
EGSnrc uses the values of the earlier 1978 edition, as EGS4 [35] did. 

MCNP, MCNPX and ITS use the electron binding energies compiled by Carlson. 

The Penelope 2008 version uses Carlson's compilation of binding energies; earlier versions 
used the compilation included in the 1978 edition of the Table of Isotopes. 
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The Geant4 toolkit uses various collections of binding energies. The main reference for bind- 
ing energies in Geant4 is the G4AtomicShells class in the materials package; according to com- 
ments in the code implementation, the binding energies values in it derive from Carlson's compi- 
lation and the 73rd edition of the CRC Handbook of Chemistry and Physics [36]. 

EADL values are used by the implementations of photon and electron interactions in Geant4 
low energy electromagnetic package [37, 38] based on the so-called Livermore Library, which 
encompasses the Evaluated Electron Data Library (EEDL) [39], the Evaluated Photon Data Library 
(EPDL97) [40] and EADL itself. EADL binding energies are also used in the calculation of proton 
ionization cross sections described in [41] and released in Geant4 9.4 in a modified version [42] 
to address the drawbacks documented in [43]. Proton ionization cross sections for Geant4 PKE 
(Particle Induced X-ray Emission) simulation described in [43] derive from ISICS [44] tabulations 
using Bearden and Burr's binding energies. 

Geant4 includes a C-i~i- reimplementation of physics models originally implemented in Pene- 
lope; Geant4 9.4 reimplements models from the 2008 version of Penelope, while previous Geant4 
releases included models equivalent to Penelope 2001. The Geant4 9.4 reimplementation appears 
to use EADL values instead of the binding energies used by Penelope 2008. Binding energies cor- 
responding to the values in the 1978 edition of the Table of Isotopes are included in a Geant4 9.4 
data set associated with Penelope. 

Ionization energies consistent with those reported by NIST [28] are included in the Geant4 
G4StaticSandiaData class. 

No reference to atomic binding energies can be retrieved in GEANT 3 documentation; how- 
ever, according to comments embedded in the code, GEANT 3 used Bearden and Burr's binding 
energies, with updated values for xenon derived from [46]. Nevertheless, the GSHLIN subroutine, 
where binding energies are hard-coded, exhibits some discrepancies with respect to both Bearden 
and Burr's tabulations and the values in [46]; the origin of these values could not be retrieved in 
the literature, nor in the software documentation. Presumably, the code implementation and its 
comments went out of phase at some stage of GEANT 3 evolution. 

Atomic binding energies are relevant to PIXE calculations; two well known software systems 
pertinent to this domain are GUPIX [47] and ISICS [44]. GUPIX uses Sevier's 1979 compilation of 
binding energies [48], which includes extensions to the 1972 collection by the same author. ISICS 
uses Bearden and Burr's binding energies by default; the most recent version of the code [49] offers 
the option of using the binding energies assembled in Williams's compilation instead of Bearden 
and Burr's ones. 

The authors of this paper could not retrieve track of the electron binding energies used by 
FLUKA in the related software documentation and in the Uterature, nor from direct inquiries with 
the maintainers of the code; it was not possible to ascertain them from the software implementation, 
whose disclosure is subject to restrictions, as their presumed source file is in a binary encoded 
format. 

2.3 Comparison of binding energies compilations 

The binding energies collected in the various compilations exhibit some differences, apart from 
those due to different references - the vacuum potential or the Fermi level. 
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A few examples of comparison are displayed in figures the plots show the difference 
between the binding energies in the various compilations and the values in Williams' compilation. 
The choice of Williams' compilation as a reference for plotting differences is arbitrary; the main 
qualitative features of the plots are anyway equivalent, if other empirical compilations are chosen 
as a reference instead of Williams' one. The difference between EADL and Williams' K shell 
binding energies is plotted separately from the other compilations, since the scale is approximately 
a factor 30 larger. 

The differences are of the order of a few electronvolts across the various empirical compila- 
tions, as illustrated in figure l(a)| , while they are larger between EADL and the empirical compila- 
tions, especially for inner shells, as shown in figure |l(b)| ; they can reach a few hundred electron- 
volts for the K shell of heavier elements. The empirical compilations derive from a common source 
(Bearden and Burr's review); therefore it is not surprising that they exhibit some similarities. 
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Figure 1. Difference between K shell binding energies in various compilations and binding energies in 
Williams' one versus atomic number: left, Carlson (blue squares). Table of Isotopes 1996 (black down 
triangles). Table of Isotopes 1978 (green up triangles), Sevier 1979 (pink stars), Bearden and Burr (empty 
circles); right, EADL (note the different scale). 



3. Strategy of the study 

The study documented in this paper is driven by pragmatic motivations. The analysis is focused on 
quantifying the accuracy of binding energy compilations used in Monte Carlo systems, and their 
impact on physics models of particle transport and on experimental observables produced by the 
simulation. The evaluation aims at identifying one or more optimal options for Monte Carlo ap- 
plications in experimental practice. A comprehensive review of the physical ground, experimental 
measurements and methods of calculations of electron binding energies is outside the scope of this 
paper. 

Two complementary approaches are adopted: direct validation of tabulated electron binding 
energy values and the evaluation of effects on related physics quantities, like ionization cross sec- 
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Figure 2. Difference between L sub-shell binding energies in various compilations and binding energies in 
Williams' one versus atomic number: EADL (red circles), Carlson (blue squares), Table of Isotopes 1996 
(black down triangles). Table of Isotopes 1978 (green up triangles), Sevier 1979 (pink stars), Bearden and 
Burr (empty circles). 



tions and X-ray fluorescence emission. Both analyses involve comparisons with experimental data 
and a comparative appraisal of the accuracy of the various compilations. 

The evaluations of effects on dependent simulated quantities represent a significant set of 
scenarios to characterize the problem domain, although they are not intended to exhaust all the 
possible present and future uses of atomic binding energies in radiation transport codes. The test 
cases are chosen to address simulation observables that depend directly and explicitly on electron 
binding energies; this requirement is essential to avoid the risk of introducing a possible bias in 
the evaluation process due to other physics modeling features than the binding energies subject to 
assessment. The test cases concern quantities that can be directly measured, so that the simulation 
outcome can be compared with experimental data to evaluate the contribution of different atomic 
binding energies to the achieved accuracy. Theoretical quantities that cannot be directly measured 
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Figure 3. Difference between representative M and N shell binding energies in various compilations 
and binding energies in Williams' compilation versus atomic number: EADL (red circles), Carlson (blue 
squares). Table of Isotopes 1996 (black down triangles). Table of Isotopes 1978 (green up triangles), Sevier 
1979 (pink stars), Bearden and Burr (empty circles). 

are explicitly not considered, since possible effects due to electron binding energies would be 
difficult to univocally identify and quantify in the analysis. 

The study is articulated over a variety of test cases, which involve different physics issues and 
reference data; the analysis methods are tailored to the physical and experimental features of each 
test case. Various statistical tools are exploited to quantify the accuracy of the distributions exam- 
ined in this study and the difference (or equivalence) of the various binding energy compilations; 

Goodness-of-fit tests mentioned in the following sections utilize the Statistical Toolkit [52,53]. 
Whenever applicable, multiple goodness-of-fit tests are applied to mitigate the risk of systematic 
effects in the conclusions of the analysis due to peculiarities of the mathematical formulation of the 
various methods. 

A combination of Student's t-test and F-test is applied to study the distribution of differences 
between the data subject to evaluation and reference values, when goodness-of-fit tests do not 
exhibit adequate discriminant power over some analyzed data samples. The t-test is utilized to 
estimate the compatibility with null mean difference, while, once the sample exhibiting the nar- 
rowest distribution of differences (i.e. the lowest variance) has been singled out, the F-test is used 
to identify the data samples whose variance is statistically equivalent to the narrowest distribution. 

The binding energies listed in the various compilations are given with respect to different 
references (vacuum or Fermi level). In the following comparisons the original values are corrected 
to account for the work function as appropriate to ensure a consistent reference. Values of the 
work function are taken from the compilation of the CRC Handbook of Physics and Chemistry 
[21], which is considered an authoritative source for these data in experimental practice; they are 
complemented by data from [54] and [55] for elements not included in the compilation of [21]. 

The analyses reported in the following sections concern elements with atomic number between 
1 and 92, unless differently specified. This range ensures uniform treatment of the various com- 
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pilations in their comparative appraisal, since all the examined compilations cover these elements, 
while only a subset of them deal with transuranic elements. Moreover, established experimental 
references of transuranic elements suitable for the analysis of binding energies are scarce. 



4. Evaluation of reference binding energies 

Comparison with experimental data is the prime method to evaluate the accuracy of simulation 
models; this validation method requires reliable experimental measurements as a reference. Three 
authoritative collections of experimental binding energies (two of which are partially overlapping) 
are used for this purpose; they include values for a limited number of elements and shells, therefore 
they can validate only part of the content of the compilations mentioned in the previous sections. 

4.1 Comparison with high precision reference data 

Experimental values of elemental binding energies reported in the literarature exhibit significant 
discrepancies [56]; they can affect the validation of binding energy compilations. 

Inconsistencies in the measurements are mostly due to inadequacies in the calibration of bind- 
ing energy scales of the various instruments, and are often visible when comparing measurements 
performed by different laboratories. 

Binding energy measurements may differ also for physical reasons: elemental binding ener- 
gies differ in the vapour and condensed states, and a chemical shift is present when atoms are 
investigated in chemical compound states. Moreover, binding energies for atoms implanted by ion 
bombardment into a metal foil substrate are shifted with respect to those for a foil of the pure el- 
ement. Measurements on different surfaces of a crystal can result in different ionization energy 
values. 

These effects may be sources of systematic errors, which can be significant when comparing 
the binding energies collected in the various compilations with experimental values. 

A further source of uncertainties derives from the conversion between binding energy values 
of solids referenced to the Fermi level and to the vacuum level; this operation involves adding, or 
subtracting, the value of the work function. 

A collection of high precision binding energies [57] was assembled by Powell for the purpose 
of constituting a reference for the NIST X-ray Photoelectron Spectroscopy Database [58]. Data 
published by different laboratories were subject to a retroactive calibration procedure; the origi- 
nal experimental values were corrected to produce a set of 61 binding energy values, concerning 
elements with atomic numbers between 4 and 84 and shells from K to N. The uncertainty of the 
reference energies is reported to be 0.061 eV [57]. These high precision data have been used to 
evaluate the accuracy of the binding energy compilations examined in this study. 

The difference of the binding energies in the various compilations with respect to the reference 
values of [57] is shown in figure 4(a)| ; the relative difference with respect to the same data is 



shown in figure 4-(b). EADL binding energies appear less accurate than the samples from other 
compilations and exhibit a systematic shift with respect to the reference values. 

It is worthwhile to remark that no quantitative, systematic validation of EADL binding ener- 
gies has been reported yet in the literature. EADL documentation [17] states in a section devoted 
to the accuracy of the data that "by comparing subshell parameters from a number of different 
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sources, it can be seen that there is still a disagreement of about 1% between the binding energies"; 
nevertheless this statement is not supported by any objective demonstration, either directly in the 
documentation itself or through references to the literature. EADL documentation does not specify 
whether the other sources considered in the above statement were other theoretical calculations, or 
empirical compilations, or experimental data. 




Figure 4. Comparison of binding energies in various compilations and reference data from [57] versus 
atomic number: EADL (red circles), Carlson (blue up triangles). Table of Isotopes 1996 (black squares). 
Table of Isotopes 1978 (green down triangles), Williams (pink stars), Sevier 1979 (turquoise asterisks). 

Goodness of fit tests, which are commonly applied in statistical analysis to compare data 
distributions, do not appear adequate to discriminate the compatibility of the various compilations 
with respect to the reference data of [57]. 

Goodness of fit tests based on the empirical distribution function (Kolmogorov-Smirnov [59, 
60], Anderson-Darling [61,62] and Cramer-von Mises [63,64]) result in p-values greater than 0.999 
for all the data samples: this means that the differences between the binding energies of the various 
compilations and the reference data are small with respect to the sensitivity of these tests to detect 
discrepancies in the distributions subject to comparison. It is worthwhile to remark that the power 
of goodness-of-fit tests is still a subject of research in statistics. 

On the other hand, the statistic [65] based on the uncertainties of the reference data reported 
in [57] (0.061 eV) results in p-values much smaller than 0.001 for all the binding energy samples; 
therefore the test would reject the hypothesis of compatibility of any compilation with the 
reference data sample with 0.001 significance. 

The outcome of the x^ test depends critically on the correct estimate of the uncertainties of the 
data subject to test. The procedure applied in [57] to build the reference data sample mitigates the 
risk of possible systematic effects due to instrumental calibration, which may affect raw experimen- 
tal data; nevertheless, other sources, independent from the intrinsic precision of the measurement, 
may contribute to the overall uncertainty. 

Previously mentioned physical and chemical shifts of the experimental data, associated with 
the conditions of the measurements, may introduce systematic effects. In this context, one should 
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take into account that, while the reference values in [57] reflect the experimental configuration 
and instrumental energy resolution for each measured material, the data tabulated in the various 
empirical compilations are the result of manipulations, such as interpolations, extrapolations and 
fits, over large collections of heterogeneous experimental data from multiple sources: the generic 
binding energy estimates deriving from these procedures may not adequately account for the pecu- 
liarities of experimental mesurements performed in specific physical and chemical configurations. 
The calculation of the test statistic includes only the uncertainties associated with Powell's 
reference data; it does not account for errors associated with the binding energies of the various 
compilations. 

A further source of uncertainty is associated with the work function in cases where a con- 
version between the Fermi and vacuum reference level should be applied for consistency between 
the distributions subject to comparison. Moreover, experimental values of the work function are 
affected by the technique of measurement and the cleanliness of the surface. The CRC compila- 
tion does not report the uncertainties of the work functions; therefore this additional error cannot 
be included in the computation of the test statistic concerning Carlson's and EADL binding 
energies. 

Due to these considerations, caution should be exercised in interpreting the outcome of the 
test as physically significant, as the nominal uncertainties values involved in the calculation may 
not realistically represent the actual uncertainties associated with the tested data samples. 

Other statistical methods than goodness-of-fit testing were exploited to quantify the compati- 
bility between the various binding energies compilations and the reference data of [57]. 

A Student's t-test was performed to estimate whether the differences between the binding 
energies of the various compilations and the corresponding reference values are compatible with a 
true mean of zero. The p- values resulting from this test aie summarized in Table [T[ The binding 
energies of Williams' and 1979 Sevier's compilations are compatible at 0.05 level of significance 
with null mean difference with respect to the reference data; the t-test rejects the hypothesis of 
compatibility with zero mean difference with 0.01 significance for all the other compilations. It 
should be stressed that these tests, as well as similar ones reported in the following, do not compare 
the compilations with each other, but how well they reproduce the set of precision reference data. 

All the binding energies compilations exhibit rather large differences with respect to the refer- 
ence data for rare gases (Ar, Xe and Kr) reported in [57], that derive from implants in other media. 
If one excludes these data from the t-test, also Bearden and Burr's binding energies are compatible 
with zero mean difference with respect to the reference data at 0.05 significance level. 

The binding energies of Williams' compilation exhibit the narrowest distribution of differences 
with respect to the reference data of [57], as can be seen in figure 4(b^ The standard deviations 
related to the various compilations are listed in Table ^, excluding the data for argon, xenon and 
krypton, which are treated as outliers. The table also reports the p-values of the F-test to evaluate 
the hypothesis of equality of variance associated with the various compilations with respect to 
the binding energies of Williams' compilation; the distributions subject to the F-test concern the 
difference between the binding energies in the compilations and the reference data of [57]. The 
statistical analysis confirms the qualitative evidence of figure 4(b)| , since the null hypothesis is 
rejected with 0.01 significance level for all the test cases. 
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Table 1. Student's t-test applied to the difference with respect to Powell's reference data. 





All data 


Excluding 


Ar, Xe, Kr 


Compilation 


Mean (eV) 


p-value 


Mean (eV) 


p-value 


Bearden and Burr 


0.51 


0.026 


0.40 


0.083 


Carlson 


0.95 


0.001 


0.65 


0.008 


EADL 


4.85 


<0.0001 


4.75 


<0.0001 


Sevier 1979 


0.41 


0.181 


0.11 


0.673 


Tol 1978 


1.00 


0.0004 


0.69 


0.002 


Tol 1996 


0.99 


0.001 


0.69 


0.006 


Williams 


0.41 


0.076 


0.09 


0.545 



Table 2. F-test applied to the differences with respect to Powell's reference data. 



Compilation 


Standard deviation (eV) 


p-value 


Bearden and Burr 


1.71 


0.0005 


Carlson 


1.80 


0.0001 


EADL 


1.56 


0.005 


Sevier 1979 


2.01 


<0.0001 


Tol 1978 


1.61 


0.002 


Tol 1996 


1.84 


0.0001 


Williams 


1.07 





4.2 Comparison with NIST recommended binding energies 

A similar analysis has been performed with respect to the collection of recommended binding 
energies for principal photoelectron lines assembled by NIST [66]. This collection consists of 85 
values; it includes most of the reference binding energies discussed in [57], along with additional 
data, mainly concerning outer shells than those reported in [57]. The data for noble gases listed 
in [57] are not included in this set of recommended values. 

Figure ^ shows the difference between the binding energies of the various compilations and 
the NIST recommended values. 

The comparison with these reference data adopts a similar method to the one described in the 
previous section. 

The p-values resulting from Student's t-test for compatibility with zero mean difference with 
respect to the reference data are reported in Table |3[ The binding energy samples extracted from the 
Williams' and Sevier's 1979 compilations are compatible with zero mean difference at 0.05 level 
of significance; those from the two editions of the Table of Isotopes are compatible at 0.01 level. 
Similarly to the previous case, the distribution of differences with the lowest standard deviation is 
the one associated with Williams' compilation. 

The standard deviations of the distributions of differences from the reference data are reported 
for all the data samples in Table ^ together with the p-values of the F-test for the equality of vari- 
ance with respect to the distribution associated with Williams' compilation. Consistently with the 
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Difference (eV) 



Figure 5. Difference between binding energies in various compilations and reference data from [66]: Carl- 
son (dashed blue line), Table of Isotopes 1996 (dotted black line). Table of Isotopes 1978 (dotted-dashed 
green line), Williams (thin pink solid line), Sevier 1979 (thick turquoise dashed line) and EADL (thick solid 
red line). 

Table 3. Student's t-test applied to the difference with respect to NIST recommended binding energies. 



Compilation 


Mean (eV) 


p-value 


Bearden and Burr 


0.35 


0.070 


Carlson 


0.75 


0.0003 


EADL 


3.96 


<0.0001 


Sevier 1979 


-0.19 


0.387 


Tol 1978 


0.44 


0.017 


Tol 1996 


0.52 


0.012 


Williams 


0.11 


0.246 



qualitative features of figure ^, the variance associated with the other compilations is incompatible 
with the variance related to Williams' binding energies. 

Based on this statistical analysis, one can conclude that Williams' compilation best reproduces 
experimental reference binding energies. It should be stressed, however, that the NIST reference 
sample represents a small subset of the periodic system of elements: approximately 6% of the total 
number of shells of elements with atomic number up to 92. 

4.3 Evaluation of ionization energies 

The ionization energy (in the past referred to as ionization potential), is the least energy that is 
necessary to remove an electron from a a free, unexcited, neutral atom, or an additional electron 
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Table 4. F-test applied to the difference with respect to NIST recommended binding energies. 



Compilation 


Standard deviation (eV) 


p-value 


Bearden and Burr 


1.74 


<0.0001 


Carlson 


1.94 


<0.0001 


EADL 


4.14 


<0.0001 


Sevier 1979 


2.03 


<0.0001 


Tol 1978 


1.66 


<0.0001 


Tol 1996 


1.86 


<0.0001 


Williams 


0.88 





from an ionized atom. In the following analysis it is considered to be equal to the binding energy 
of the least bound electron in the atom. 

A compilation of reference experimental ionization potentials is available from NIST [28] ; the 
same values are also reported in the Table of Isotopes [18] in a table distinct from the compila- 
tion of electron binding energies and in the CRC Handbook of Chemistry and Physics [21]. This 
compilation does not list the uncertainties of the ionization energies it collects, but NIST web site 
comments that they range from less than one unit in the last digit of the given values to more than 
0.2 eV. 

The lowest binding energies for each element in the various compilations have been compared 
to the reference ionization energies collected by NIST. The compilations of the 1978 edition of the 
Table of Isotopes and Williams do not include many outer-shell binding energies; this limitation 
may be related to the emphasis of these compilations on experimental effects related to inner shells, 
like measurements concerning X-ray fluorescence or Auger electron emission. Therefore the fol- 
lowing analysis was restricted to the compilations of Bearden and Burr, Carlson, EADL, Sevier 
1979 and 1996 edition of the Table of Isotopes, which list a full set of electron binding energies. 

The difference between ionization energies derived from the various compilations and NIST 
reference data is shown in figure ^. Carlson's binding energies appear to be in closest agreement 
with NIST ionization energies. 

Goodness of fit tests are sensitive to the differences exhibited by the various compilations 
with respect to the NIST reference collection. The results of the Kolmogorov-Smimov, Anderson- 
Darling and Cramer-von Mises tests are listed in Table ||. The hypothesis of compatibility with 
NIST reference data is rejected by all the tests with 0.001 significance for EADL and Bearden 
and Burr's data. Carlson's compilations and the Table of Isotopes 1996 are compatible with the 
reference data at 0.05 significance level according to all the tests. Regarding Sevier's 1979 compi- 
lation, the Anderson-Darling test rejects the hypothesis of compatibility at 0.05 significance level, 
while the Kolmogorov-Smimov and Cramer-von Mises test do not; the different response of these 
tests near the critical region of 0.05 significance could be explained by the greater sensitivity of the 
Anderson-Darling test statistic to fat tails. 
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Difference (eV) 

Figure 6. Difference between ionization energies in various compilations and NIST reference experimental 
data: Carlson (thin solid blue line), Table of Isotopes 1996 (dashed black line), Sevier 1979 (thick solid 
turquoise line) and EADL (dot-dashed red line). 

Table 5. P-values from goodness-of-fit tests concerning NIST reference ionization potentials. 



Compilation 


Kolmogorov 


Anderson 


Cramer 




Smirnov 


Darling 


von Mises 


Bearden and Burr 


<0.001 


<0.001 


<0.001 


Carlson 


0.670 


0.995 


0.963 


EADL 


<0.001 


<0.001 


<0.001 


Sevier 1979 


0.061 


0.023 


0.067 


Tol 1996 


0.096 


0.099 


0.213 



5. Effects on fluorescence X-ray energies 

Compilations of characteristic X-ray energies are available, at least for lines of experimental in- 
terest, which in principle could be used in Monte Carlo simulation to determine the energy of 
secondary products of atomic relaxation. Nevertheless, these experimental tabulations can hardly 
satisfy the requirements of general purpose Monte Carlo codes, which require the ability of gener- 
ating any atomic transition for any element. The energies of X-rays and Auger electrons resulting 
from atomic relaxation are often computed by Monte Carlo codes as the difference between the 
binding energies of the shells involved in the transition; in this approximation the binding energies 
of the atom in an ionised state are assumed to be the same as in the ground state. Therefore the 
accuracy of the simulation of the secondary products of atomic relaxation is determined by the 
accuracy of the binding energies implemented in the Monte Carlo system (apart from the physical 
approximation of neglecting the difference between the binding energies of an ionised atom and a 



-15- 



neutral one in the ground state). 

The accuracy of the examined binding energy compilations to reproduce the energy of atomic 
relaxation products has been estimated with respect to the experimental X-ray energies reported in 
the review by Deslattes et al. [67], which concerns K and L transitions. 

A comparison of X-ray energies calculated by Geant4, based on EADL binding energies, with 
respect to the same experimental data is documented in [68]. That study showed that, according 
to the outcome of the Kolmogorov-Smirnov test, all the X-ray energies simulated by Geant4 are 
compatible with the experimental data with 0.1 significance for all transitions and all elements; 
the relative difference between simulated and experimental values is approximately 1-2% for most 
individual transitions. The present study is extended to binding energy compilations other than 
EADL. 

A selection of representative plots of the relative difference between calculated X-ray energies 
and the experimental data of [67] is shown in figures X-ray energies are calculated from the 
various compilations of binding energies. It is evident from the plots that the energies calculated 
by EADL appear less accurate than those based on the other compilations. Nevertheless, as already 
found in [68], the discrepancies of the energies deriving from EADL with respect to measurements 
are quite small (less than 2% in general). 

Similarly to what has been discussed in the previous section, goodness-of-fit tests based on 
the empirical distribution function are not sensitive to such small differences: the hypothesis of 
compatibility between experimental data and X-ray energies based on EADL (the compilation that 
is evidently responsible for the largest discrepancies) is not rejected at 0. 1 level of significance [68]. 

The test has limited discriminant power as well, due to the small uncertainties of the exper- 
imental reference data in [67] (less than 0.1 eV for some transitions), which lead to the rejection of 
the null hypothesis of compatibility between calculated and experimental X-ray energies in a large 
number of test cases. It is hard to ascertain whether this result of the test is due to underesti- 
mated uncertainties for some transitions and elements, or reflects a realistic conclusion that X-ray 
energies calculated from binding energy differences do not achieve the same accuracy by which 
X-ray energies are experimentally measured. 

Similarly to what was described in the previous section, a t-test was applied to evaluate whether 
the distribution of differences between calculated and experimental X-ray energies is compatible 
with a true mean of zero. For each transition the t-test was performed over all the elements for 
which experimental values are reported in [67]; the fraction of test cases for each transition for 
which the hypothesis of compatibility is not rejected with 0.05 significance is listed in Table ^. 

The largest number of test cases where the hypothesis of compatibility with null average dif- 
ference is not rejected with 0.05 significance is achieved by the compilations of the 1996 Table of 
Isotopes and Bearden and Burr's review (44 out of 48 test cases). 

The hypothesis whether the compatibility of the other binding energy compilations with zero 
mean is equivalent to the one achieved by the 1996 Table of Isotopes and Bearden and Burr's 
review was tested by means of contingency tables. Contingency tables were built by counting in 
how many t-test cases the rejection of the null hypothesis occurs, or does not occur; these counts 
are respectively identified as "fail" or "pass". The results concerning K and L shells are summed to 
obtain a larger sample size. They were analyzed by means of Fisher's exact test [70], Pearson's 
test (whenever the number of entries in each cell justifies the use of this test) and the test with 
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Table 6. Number of test cases compatible at 0.05 significance level with mean null difference between 
calculated and experimental X-ray energies. 



Compilation 


K 


Li 


L2 


L3 


K+L 


Fraction 


Bearden and Burr 


10 


12 


11 


11 


44 


0.92 ± 0.04 


Carlson 


7 


11 


11 


9 


38 


0.79 ± 0.06 


EADL 


8 


4 


1 


3 


16 


0.33 ± 0.07 


Sevier 1979 


10 


9 


12 


9 


40 


0.83 ± 0.05 


Tol 1978 


7 


11 


8 


9 


35 


0.76 ± 0.06 


Tol 1996 


9 


12 


11 


12 


44 


0.92 ± 0.04 


Williams 


8 


12 


11 


9 


40 


0.85 ± 0.05 



Table 7. Contingency table associated with the t-test : applied to X-ray energies derived from EADL and 
from the 1996 Table of Isotopes. 



'X}^ test outcome 


Tol 1996 


EADL 


Pass 


44 


16 


Fail 


4 


32 


p-value Fisher test 


< 0.0001 


p-value Pearson 'x} 


not applicable 


p-value Yates %^ 


< 0.0001 



Yates continuity correction [72]. The contingency table concerning the comparison of EADL and 
the 1996 Table of Isotopes is reported in Table 0, along with the p- values of the three tests applied 
to it. 

The hypothesis of equivalence with respect to the results of the t-test is rejected with 0.05 
significance for EADL; it is not rejected for Carlson's, Sevier's and Williams' compilations. The 
outcome of the tests is controversial for the contingency table concerning the 1996 and 1978 edi- 
tions of the Table of Isotopes: the p-values are 0.050 for Fisher's exact test, 0.039 for Person's 
test and 0.075 for the with Yates' continuity correction. The compatibility between EADL and 
the the 1996 Table of Isotopes is excluded even if a looser 0.01 significance for the rejection of the 
null hypothesis is set both in the t-test and in the contingency tables. 

The distribution of the difference between the X-ray energies calculated from binding energy 
tabulations and the experimental values of [67] is wider for EADL than for all the other compila- 
tions; this result can be appreciated in a few representative plots (figures |-12). 

The equivalence of the variance of the differences between calculated and experimental X-ray 
energies was estimated by means of the F-test. For each transition, the variance of the correspond- 
ing data sample was compared to the variance associated with the 1996 Table of Isotopes, which 
exhibits the narrowest distribution of differences between calculated and experimental X-ray en- 
ergies. The hypothesis of equivalence of the variances under test was rejected if the p-value from 
the F-test was smaller than 0.01. The fraction of transitions for which the outcome of the F-test 
indicates that there is no significant difference in the respective variances is listed in Table ||. 
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Figure 7. K-shell transitions, relative difference between binding energies in various compilations and 
experimental data from [67] versus atomic number: EADL (red circles), Carlson (blue up triangles). Table 
of Isotopes 1996 (black squares). Table of Isotopes 1978 (green down triangles), Williams (pink stars), 
Sevier 1979 (turquoise asterisks) and G4AtomicShells (empty squares). 

Table 8. Fraction of K and L transitions for which the variance of the difference between calculated and 
experimental X-ray energies is equivalent to the variance associated with the 1996 Table of Isotopes. 



Compilation 


Fraction of transitions 


Beaiden and Burr 


0.60 ± 0.07 


Carlson 


0.69 ± 0.07 


EADL 


0.25 ± 0.06 


Sevier 1979 


0.79 ± 0.06 


Tol 1978 


0.65 ± 0.07 


Williams 


0.70 ± 0.07 
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Figure 8. Li -shell transitions, relative difference between binding energies in various compilations and 
experimental data from [67] versus atomic number: EADL (red circles), Carlson (blue up triangles). Table 
of Isotopes 1996 (black squares). Table of Isotopes 1978 (green down triangles), Williams (pink stars), 
Sevier 1979 (turquoise asterisks) and G4AtomicShells (empty squares). 
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Figure 9. L2-shell transitions, relative difference between binding energies in various compilations and 
experimental data from [67] versus atomic number: EADL (red circles), Carlson (blue up triangles). Table 
of Isotopes 1996 (black squares). Table of Isotopes 1978 (green down triangles), Williams (pink stars), 
Sevier 1979 (turquoise asterisks) and G4AtomicShells (empty squares). 



The results of the F-test are consistent with the qualitative appraisal of the accuracy of the 
distributions in figures |5 - 12. It is worthwhile to recall that the F-test is sensitive to the normality of 
the distributions to which is applied; although the differences between calculated and experimental 
data are expected to be normally distributed, the results reported in Table ^ may be affected by 
some details of the distributions subject to comparison. 

The analysis of X-ray energies suggests that better accuracy in the reproduction of K and L 
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Figure 10. L3-sheU transitions, relative difference between binding energies in various compilations and 
experimental data from [67] versus atomic number: EADL (red circles), Carlson (blue up triangles). Table 
of Isotopes 1996 (black squares). Table of Isotopes 1978 (green down triangles), Williams (pink stars), 
Sevier 1979 (turquoise asterisks) and G4AtomicShells (empty squares). 
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Figure 11. -shell transitions, relative difference between binding energies in various compilations and ex- 
perimental data from [67]: EADL (thick solid red line), Carlson (dashed blue line). Table of Isotopes 1996 
(thin solid black line), Williams (dash-dotted green line), Sevier 1979 (dotted turquoise line); the results of 
the other compilations considered in this study, which are not shown, exhibit a narrow distribution similar to 
the other compilations, except EADL. 



transition energies can be achieved by binding energy compilations other than EADL. A comment 
on the accuracy of EADL transition energies of the M series, mentioning deviations from experi- 
ment up to 10 MeV, is reported in the X-Ray Data Handbook [9] with the recommendation of using 
other experimental data preferably; nevertheless, no comprehensive, quantitative demonstration of 
EADL accuracy is reported, nor references are cited in support of the appraisal of EADL accuracy 
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Figure 12. L-shell transitions, relative difference between binding energies in various compilations and 
experimental data from [67]: EADL (thick solid red line), Carlson (dashed blue line). Table of Isotopes 
1996 (thin solid black line), Williams (dash-dotted green line), Sevier 1979 (dotted turquoise line); the 
results of the other compilations considered in this study, which are not shown, exhibit a narrow distribution 
similar to the other compilations, except EADL. 



and consequent recommendation. 

The original design of Geant4 atomic relaxation described in [69] would easily accommodate 
the improvement of the accuracy of the simulated energies through alternative binding energy op- 
tions: the software implementation would handle the process of atomic relaxation transparently, if 
a different tabulation of binding energies is supplied as an external file. 



6. Effects on ionization cross sections 



Some analytical formulations of cross sections for the ionization of atoms by charged particle im- 
pact involve atomic binding energies. Two of these models are considered in this study to ascertain 
whether different binding energy compilations would produce significant differences in the cross 
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section values: the Binary-Encounter-Bethe model (BEB) [73] for electron impact ionization and 
the ECPSSR (Energy Loss Coulomb Repulsion Perturbed Stationary State Relativistic) model [74] 
for proton impact ionization. For both models the effects on the accuracy of the cross section 
calculations are quantitatively estimated through a comparison with experimental data. 

6.1 Electron impact ionization cross sections 

Two models of electron impact cross sections, the Binary-Encounter-Bethe [73] model and the 
Deutsch-Mark model [75], have been designed, implemented and validated in view of extending 
and improving Geant4 simulation capabilities in the energy range below 1 keV. Their features, 
verification and validation are briefly summarized in [76, 77] and extensively documented in a 
dedicated paper [78]. The first software development cycle has been focussed on modeling total 
ionization cross sections; the validation process and the analysis of the effect of atomic binding 
energies concern these calculations, although the BEB model has the capability of calculating cross 
sections for the ionisation of individual shells. 

The BEB cross section for the ionization of subshell / is given by: 



t + {u+\)/n 
where: 



T U ^ . ,.JR 



(6.1) 



t = -, u = -, S = 4nat,N\^-j (6.2) 

In the above equations B is the electron binding energy, is the the occupation number, T is the 
incident electron energy, U is the average electron kinetic energy, t and u are normalized incident 
and kinetic energies, n is the principal quantum number (only taken into account when larger than 
2), ao is the Bohr radius and R is the Rydberg constant. The sum over all the subshells / of an atom 
gives the total (counting) cross section; in practice, only the valence shell and a few outer subshells 
contribute significantly to determine the cross section value. 

The original BEB cross sections [73] used binding energy values calculated by the authors of 
the model. Only a few of those values are documented in [73] ; they were utilized in the software 
verification process to assess the correctness of the implementation, but such a small set is inade- 
quate for using the model in a general purpose simulation system, which must be able to calculate 
cross sections for any target atoms. The BEB model developed for use with Geant4 utilizes a full 
set of binding energies and provides the option of accessing alternative compilations. 

The analysis addressed two issues: the sensitivity of cross sections to the values of the binding 
energies used in the calculation, and the evaluation of the accuracy with respect to experimental 
data. 

Two examples of the effects of different binding energies on the calculated cross sections are 
shown in figure |13[ They illustrate three options of binding energies: Lotz's compilation, which is 
also used by the Deutsch-Mark model, EADL data for all shells and EADL with ionization poten- 
tials replaced by NIST values [28] (identified in the following as "modified EADL"). Lotz's com- 
pilation is identical to Carlson's apart from a few exceptions; according to the analysis in section 



4.3, Carlson's compilation appears the most accurate with respect to NIST ionization potentials. 



while EADL exhibits the largest differences with respect to other compilations in both inner and 
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outer-shell binding energies. Significant differences are visible in the cross sections, when different 
ionization energies are used in the calculation, while different inner-shell binding energies appear 
to have relatively small effects. 



£ 2-5 



Energy (eV) 

(a) Nitrogen 




Energy (eV) 

(b) Germanium 



Figure 13. BEB electron impact ionization cross section, BEB model with EADL binding energies for all 
shells (empty circles), BEB model with EADL binding energies except for ionization energies replaced by 
NIST values (empty squares), BEB model with Lotz binding energies (asterisks) and experimental data: (a) 
from [106] (red squares), [107] (black circles) and [108] (green stars), (b) from [121] (black circles). 

The effect of these three options, which can be considered as extreme alternatives in the BEB 
calculation, has been quantified through a statistical analysis. 

First, cross sections calculated with different binding energies were compared via goodness- 
of-fit tests; the test concerned all elements with atomic number between 1 and 92, and incident 
electron energies from 1 eV to 10 keV, divided in two ranges: those up to 100 eV, and those above. 
The results of the comparisons are summarized in Tables^ and 10. The Kolmogorov-Smirnov test 
appears the most sensitive to differences in the cross sections deriving from the considered binding 
energy options. The hypothesis of compatibility between cross sections calculated with Lotz's 
binding energies and with modified EADL is rejected with 0.05 significance only for a few heavy 
elements (with Z>80) in the lower energy range (below 100 eV): this result indicates that total 
ionization cross sections are marginally affected by inner-shell binding energies. The hypothesis 
of compatibility between the cross sections based on EADL and modified EADL is rejected in a 
larger number of test cases, especially in the lower energy range: this result shows that the cross 
sections are sensitive to the values of the ionization potential. 

The following analysis evaluated whether different ionization potentials would significantly 
affect the accuracy of the calculated cross sections with respect to experimental data. The effects of 
different ionization energies on the accuracy of the calculation are not straightforward to ascertain 
from a qualitative appraisal of the data: in fact, within the data sample one can identify test cases 
where either configuration - with NIST values or with EADL original values - appears to better re- 



produce the experimental data as shown, for instance, in figure [13|. Therefore a statistical analysis 
was performed, examining the compatibility with experiment of two sets of BEB cross sections. 
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Table 9. Test cases in which the hypothesis of compatibility of BEB cross sections based on Lotz and 
modified EADL binding energies is not rejected. 



Test 


Fraction 

(E < 100 eV) 


Fraction 

(E > 100 eV) 


Kolmogorov-Smirnov 
Anderson-Darling 
Cramer-von Mises 


0.91 ± 0.05 
0.97 ± 0.03 
0.97 ± 0.03 


1.00-0.02 
1.00-0.02 
1.00-0.02 



Table 10. Test cases in which the hypothesis of compatibility of BEB cross sections based on EADL and 
modified EADL is not rejected. 



Test 


Fraction 

(E < 100 eV) 


Fraction 

(E > 100 eV) 


Kolmogorov-Smirnov 
Anderson-Darling 
Cramer-von Mises 


0.70 ± 0.05 
0.83 ± 0.04 
0.95 ± 0.02 


0.98 ± 0.02 
0.98 ± 0.02 
1.00 - 0.02 



which use respectively EADL binding energies for all shells, or the modified EADL with NIST 
ionization potentials. The two sets of cross sections were compared to the same experimental mea- 
surements [79]- [168], consisting of more than 120 individual data sets and concerning more than 
50 elements. The comparison with experimental data exploits goodness-of-fit tests (Kolmogorov- 
Smirnov, Anderson-Darling, Cramer-von Mises and X^); their significance was set to 0.05. The 
test was articulated over five distinct energy ranges below 1 keV to appraise in detail the accuracy 
of the calculated cross sections. 

The number of test cases for which the null hypothesis of compatibility between calculated and 



measured cross sections is rejected, or not rejected, is reported in Table 1 1 for the two examined 
binding energy options. Below 1 keV the hypothesis of compatibility with experimental data is 
always rejected in a smaller number of test cases when the cross section calculation utilizes NIST 
ionization energies instead of EADL original ones. No difference is observed above 1 keV. 

The analysis by means of contingency tables does not reject the hypothesis of equivalence 
between the two cross section categories at reproducing experimental data in any of the considered 
energy ranges. Nevertheless, the probability that the better performance associated with NIST 
ionization energies in all five trials could be due to chance only is 0.03. 

6.2 Proton impact ionization cross sections 

A similar study was performed on proton ionization cross sections. Several cross section models 
for the computation of inner-shell ionization by proton and a particle impact have been released 
in Geant4 version 9.4 [43]; they include calculations based on the plane wave Bom approximation 
(PWBA) [169], the ECPSSR (Energy-loss Coulomb Perturbed Stationary State Relativistic) model 
[74] in a number of variants and a collection of empirical models, deriving from fits to experimental 
data. The PWBA and ECPSSR cross sections (in all their variants) exploit tabulations produced by 
the ISICS (Inner-Shell Ionization Cross Sections) code [44] for K, L and M shells. 
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Table 11. Test cases for which the hypothesis of compatibility of BEB cross sections with experimental data 
is rejected, or not rejected, with 0.05 significance. 



Energy 




EADL 




Modified EADL 


(eV) 


Pass 


Fail 


Pass Fraction 


Pass 


Fail 


Pass Fraction 


<20 


67 


40 


0.63±0.05 


79 


28 


0.74±0.04 


20-50 


61 


68 


0.47±0.04 


81 


48 


0.63±0.04 


50-100 


40 


84 


0.32±0.04 


49 


75 


0.40±0.04 


100-250 


47 


80 


0.37±0.04 


56 


71 


0.44±0.04 


250-1000 


45 


31 


0.59±0.06 


47 


29 


0.62±0.06 


>1000 


14 


11 


0.56±0.10 


14 


11 


0.56±0.10 



The formulation of the PWBA and ECPSSR cross sections involves atomic binding energies. 
For a given shell the PWBA cross section is given by 

cTpwBA = aoe-^Ff^,e') (6.3) 



,0 

where: 



Z 



a,, = SKat, [ ^ ) (6.4) 



ao is the Bohr radius, Zi is the projectile atomic number, Z2 is the effective atomic number of the 
target atom, F is the reduced universal cross section, with the reduced atomic electron binding 
energy 6 and reduced projectile energy rj given by 



and 



e=2«2g (6.5) 



respectively. In equations 6.5 and 6.6 E, M and U represent the energy, mass and atomic binding 
energy. In the above formulae the indices 1 and 2 refer respectively to the projectile and the 
target. The analytical formulation of the reduced universal cross section F can be found in [44] ; it 
involves the reduced atomic electron binding energy. The ECPSSR cross section for a given shell 
is expressed in terms of the PWBA value: 



aECPSSR = C| {dqa^Qovm | ^^^y | (6.7) 

where C| is the Coulomb deflection correction, is the correction factor for binding energy and 
polarization effects, niR is the relativistic correction, is the minimum momentum transfer and 

^=vi^ (6.8) 

vi being the projectile velocity. 
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The PWBA and ECPSSR cross section tabulations distributed with the Geant4 code were pro- 
duced with the ISICS 2008 version, which uses the Bearden and Burr binding energies. Recent 
updates to ISICS [49] offer the option of using Williams' compilation of binding energies as alter- 
native values to the default Bearden and Burr's ones; a further evolution of ISICS [170] lets the user 
specify an arbitrary source of atomic binding energies, thus providing access to any of the options 
analyzed in this paper. This version, which was used to produce the data for this paper, involves 
new implementations of some parts of the ISICS code, which contribute to the numerical correct- 
ness and computational robustness of the software. The new features and verification of this new 
version of ISICS are documented in [170]. The experimental validation of cross sections generated 
by this new version of ISICS produces consistent results with those reported in [43], when the two 
code versions are run in the same configuration. 

Cross sections calculated by ISICS 2011 version using different binding energy compilations 
exhibit differences; some examples, concerning K and L shells, are shown in figures |I4-16. The 
differences appear larger for light elements and K shell. 




Figure 14. K shell ionization cross sections by proton impact on carbon (a) and copper (b) calculated by 
the ECPSSR model with different binding energies: the plot shows the relative difference with respect to the 
values calculated with Bearden and Burr's binding energies used by default by ISICS; the symbols identify 
cross sections calculated with EADL (red circles), Carlson (blue squares), 1996 Table of Isotopes (black up 
triangles), 1978 Table of Isotopes (green down triangles), Williams (pink empty circles) and Sevier 1979 
(turquoise stars) binding energies. Some symbols are not visible in the plot due to the close values of some 
cross sections. 

The effects of different binding energies on the accuracy of proton ionization cross sections 
have been evaluated by comparing values based on various binding energy collections with ex- 
perimental data for K and L shells. The experimental data derive from the reviews by Paul and 
Sacher [171], Sokhi and Crumpton [172] and Orlic et al. [173]; the comparison process adopts the 
same strategy described in [43] for the validation of the cross section models available in Geant4. 
For each element, the compatibility between calculated and experimental cross sections is evalu- 
ated by means of the 'y} test; the significance of the test for the rejection of the null hypothesis of 
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Figure 15. L sub-shell ionization cross sections by proton impact on cadmium calculated by the ECPSSR 
model with different binding energies: the plot shows the relative difference with respect to the values calcu- 
lated with Bearden and Burr's binding energies used by default by ISICS; the symbols identify cross sections 
calculated with EADL (red circles), Carlson (blue squares), 1996 Table of Isotopes (black up triangles), 1978 
Table of Isotopes (green down triangles), Williams (pink empty circles) and Sevier 1979 (turquoise stars) 
binding energies. Some symbols are not visible in the plot due to the close values of some cross sections. 



equivalence of the compared distributions is set to 0.05. 

The analysis of the sensitivity to electron binding energies is reported here for plain ECPSSR 
cross sections. The fraction of tested elements for which K shell cross sections calculated with 
various binding energies are compatible with experimental data is listed in Table [l^: with the ex- 
ception of EADL, all binding energy compilations appear to produce equivalently accurate cross 
sections. The use of EADL binding energies results in fewer test cases that are compatible with 
measurements; the hypothesis of equivalent accuracy of cross sections based on EADL with those 
based on other compilations is rejected with 0.05 significance for all the alternative binding ener- 
gies. As an example, the contingency table comparing the compatibility with experimental data of 
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Figure 16. L sub-shell ionization cross sections by proton impact on tungsten calculated by the ECPSSR 
model with different binding energies: the plot shows the relative difference with respect to the values calcu- 
lated with Bearden and Burr's binding energies used by default by ISICS; the symbols identify cross sections 
calculated with EADL (red circles), Carlson (blue squares), 1996 Table of Isotopes (black up triangles), 1978 
Table of Isotopes (green down triangles), Williams (pink empty circles) and Sevier 1979 (turquoise stars) 
binding energies. Some symbols are not visible in the plot due to the close values of some cross sections. 



cross sections based on EADL and Bearden and Burr binding energies is reported in Table 13 



The differences of cross sections associated with binding enegy compilations are smaller for L 
shell than for K shell, as one can qualitatively observe in two examples, concerning cadmium and 
tungsten, shown in figures 15-1^. The comparison of L shell cross sections with experimental data 
does not identify any significant differences associated with the use of different binding energies; 
the fraction of tested elements for which L shell cross sections calculated with various binding 
energies are compatible with experimental data is listed in Table 14. 

The scarcity of experimental measurements prevents a similar analysis on the effect of binding 
energies on M shell ionization cross sections. 
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Table 12. Fraction of tested elements for which ECPSSR K shell cross sections are compatible with experi- 
mental data. 



^UmpilallUIl 


P raCllUIl 


Bearden and Bun- 


0.76 ± 0.05 


Carlson 


0.76 ± 0.05 


EADL 


0.58 ± 0.06 


Sevier 1979 


0.76 ± 0.05 


Tol 1978 


0.78 ± 0.05 


Tol 1996 


0.76 ± 0.05 


Williams 


0.78 ± 0.05 



Table 13. Contingency table to estimate the equivalent accuracy of ECPSSR K shell cross sections using 
EADL and Bearden and Burr's binding energies. 



test outcome 


Bearden and Burr 


EADL 


Pass 


51 


39 


Fail 


16 


28 


p-value Fisher test 


0.042 




p-value Pearson 


0.027 




p-value Yates x^ 


0.043 





Table 14. Fraction of tested elements for which ECPSSR L shell cross sections are compatible with experi- 
mental data. 



Compilation 


Fraction, Li 


Fraction, L2 


Fraction, L3 


Bearden and Bun- 


0.68 ± 0.09 


0.68 ± 0.09 


0.89 ± 0.06 


Carlson 


0.68 ± 0.09 


0.68 ± 0.09 


0.86 ± 0.07 


EADL 


0.64 ± 0.09 


0.64 ± 0.09 


0.89 ± 0.06 


Sevier 1979 


0.68 ± 0.09 


0.68 ± 0.09 


0.89 ± 0.06 


Tol 1978 


0.68 ± 0.09 


0.68 ± 0.09 


0.89 ± 0.06 


Tol 1996 


0.68 ± 0.09 


0.68 ± 0.09 


0.89 ± 0.06 


Williams 


0.68 ± 0.09 


0.68 ± 0.09 


0.89 ± 0.06 



Based on this analysis, the accuracy of proton ionization cross sections appears statistically 
equivalent for all binding energy options but EADL. 

7. Effects on Compton scattering 

Doppler broadening of photon energy spectra arises from Compton scattering between photons 
and moving electrons bound to atoms of the target medium. Algorithms to account for Doppler 
broadening are implemented in widely used Monte Carlo systems: those included in EGS [174], 
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Figure 17. Energy distribution of photons between 89° and 91° resulting from Compton scattering of 40 
keV photons orthogonally impinging onto a silicon target, obtained using EADL (red) and Carlson's (black) 
binding energies in the simulation. The two histograms associated with either binding energy options are 
practically undistinguishable. 



MCNP [175] and Geant4 [176] are based on the method described in [174]; the algorithm imple- 
mented in Penelope produces equivalent results [176]. 

A test was performed to ascertain if different binding energy compilations would affect the cal- 
culated energy distributions of Compton scattering generated in the simulation. The test concerned 
a few target materials relevant to Compton telescopes [177], silicon, germanium and xenon, which 
are characterized by different experimental resolutions related to the effects of Doppler broadening. 
For this investigation, the original implementation of Compton scattering with Doppler broadening 
in Geant4 and associated unit test [176] was used. The analysis compared the spectra deriving from 
two sets of binding energies: those used in the simulation, which derive from EADL, and Carlson's 
compilation. The latter was chosen as its binding energies for the considered elements exhibit the 
largest difference with respect to EADL ones among the various examined compilations. 

No significant effect was visible in the spectra of the scattered photons as a result of simula- 



tions using different binding energy compilations. An example is illustrated in figure 17, which 
shows the energy spectrum of photons between 89° and 91° resulting from Compton scattering of 
40 keV photons orthogonally impinging onto a silicon target. 

Pearson's test confirms the equivalence of the Doppler broadened photon spectra based on 
EADL and Carlson's binding energies; the p-value resulting from this test is 1 for all the three 
target materials. 

Therefore, based on this investigation, one can conclude that the choice of binding energy 
compilation is not critical for the simulation of Compton scattering accounting for Doppler broad- 
ening. 

8. Merged compilations 

The collection of binding energies in Geant4's G4AtomicShells class has been assembled specifi- 



-30- 



10 20 30 40 50 60 70 SO 90 100 



Atomic number 

Figure 18. KL3 transition, difference between X-ray energies calculated from binding energies and exper- 
imental data from [67] versus atomic number: binding energies from G4AtomicShells (red circles), from 
Carlson (blue squares) and Williams (black triangles). 



cally for Geant4, merging data from Carlson's compilation with others from the 73rd edition of the 
CRC Handbook of Chemistry and Physics [36]. The origin of the data for each element and shell 
is documented in the form of comments in the code implementation. The authors of this paper 
could not retrieve a copy of the latter reference, which has been superseded by more recent editions 
(the most recent one at the time of writing this paper is the 91st edition); nevertheless, most of 
the values identified in the comments to the code as originating from [36] appear consistent with 
those published in the most recent version of the Handbook, which includes Williams' compila- 
tion. A few values in G4AtomicShells, however, are consistent with neither Carlson's nor WilUams' 
compilations. 

The two sources, Carlson's and Williams' compilations, report binding energies based on dif- 
ferent reference levels: the vacuum level for Carlson's data and the Fermi level for Williams' 
data. Data referring to different reference levels are associated with shells of the same element 
in G4AtomicShells. The inconsistency of the data in the G4AtomicShells class may generate sys- 
tematic effects in physics observables; some examples are illustrated in figures 18 and 19. These 
plots show the differences between X-ray energies calculated from G4AtomicShells binding ener- 
gies and the experimental data of Deslattes et al. [67], along with X-ray energies calculated from 
Carlson's and Williams' compilations: the X-ray energies based on G4AtomicShells exhibit some 
systematic shifts with respect to the experimental data, while the X-ray energies based on Carl- 
son's and Williams' compilations do not appear affected by such systematic discrepancies with 
measurements. The systematic effect is so large, that statistical tests appear redundant to identify 
its occurrence. 



9. Conclusion 



A survey of compilations of atomic binding energies compilations used by general purpose Monte 
Carlo transport codes and other specialized software systems has been performed. Most compila- 
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Figure 19. KM2 transition, difference between X-ray energies calculated from binding energies and exper- 
imental data from [67]: binding energies from G4AtomicShells (red shaded histogram), Carlson (blue solid 
line histogram) and Williams (black dashed line histogram). 



tions are based on experimental data; the only exception among those considered in this study is 
EADL, which is the result of theoretical calculations. 

The accuracy of these compilation has been evaluated through direct comparisons with exper- 
imental data and through their effects on related physical quantities used in particle transport and 
experimental observables. 

The results of this study show that no single compilation is ideal for all applications. 

Direct comparisons with reference experimental data, which concern a subset of K, L, M and 
N shells, identify Williams' compilation, included in the X-ray Data Booklet and CRC Handbook 
of Physics and Chemistry, as the one best agreeing with experimental values. Regarding outer 
shells, ionization energies appear to be best reproduced by Carlson's compilation; Lotz's ionization 
energies are identical to Carlson's for elements with atomic number up to 92. 

K and L shell X-ray energies are more accurately calculated based on the binding energies 
reported in the 1996 Table of Isotopes. With respect to these binding energies, the X-ray energies 
derived from the compilations by Bearden and Burr, Sevier and Williams do not exhibit any sta- 
tistically significant disagreement in compatibilility with average zero difference from experiment. 
Regarding the variance of such differences, the various compilations, with the exception of EADL, 
are equivalent to the 1996 Table of Isotopes in 60% to 79% of the transitions. X-ray energies based 
on EADL, although less accurate than those produced by other compilations, differ from the ex- 
perimental references by less than 2% for most transitions: such an inaccuracy can be tolerable in 
some experimental applications, while others, where accuracy of simulated X-ray energies is im- 
portant, should utilize the 1996 Table of Isotopes or other compilations providing better accuracy 
than EADL in this domain. 

Total cross sections for electron impact ionization addressing the energy range below 1 keV 
(relevant to microdosimetry applications) are sensitive to the values of ionization potentials, while 
they are marginally influenced by inner shell binding energies. Lotz-Carlson's compilation and 
EADL modified to include NIST experimental ionization energies exhibit equivalent behaviour. 



-32- 



while the use of EADL original ionization energies substantially decreases the accuracy of the 
cross section calculation. 

K shell ionization cross sections by proton impact are sensitive to the binding energy values 
used in the ECPSSR calculation. All the empirical binding energy compilations produce results 
compatible with experimental measurements, while cross sections based on EADL show statisti- 
cally significant worse accuracy. Binding energies also affect the calculation of ECPSSR cross 
sections for the L shell, but the differences are smaller than for the K shell; the uncertainties of L 
shell experimental data are too large to appreciate the effect of different binding energies in terms 
of modehng accuracy. 

The simulation of Doppler broadening in Compton scattering appears insensitive to the choice 
of binding energies among the examined compilations. 

It is worthwhile to stress that these quantitative conclusions were not previously documented 
in the literature. 

Simulation appHcations with high precision requirements may profit from the results docu- 
mented in this paper to identify the optimal set of binding energies for specific scenarios (e.g. 
material analysis, microdosimetry, FIXE etc.). Simulations appUcations not characterized by high 
precision requirements may be satisfied by any of the compilations. Monte Carlo kernel developers 
may profit form the results documented in this paper to optimize the accuracy of the physics models 
they deliver to the experimental user community for applications. 

As no single compilation is suitable for all applications, it is highly desirable for simulation 
packages to allow experimentalists to choose which compilation to use for their application. This 
is much more easily possible in packages which read their binding energy data from file at run- 
time, rather than hard-coding it into the appUcation binaries as in GEANT 3 and Geant4 materials 
package. 

EADL binding energies appear consistently associated with worse accuracy in all the test 
cases analysed in this paper. An evolution of EADL to better reflect the state-of-the-art would be 
desirable; it has already been advocated in [178] regarding the improvement of radiative transition 
probabilities. However, modifications to EADL should preserve the consistency with two other 
related data libraries, EEDL (Evaluated Electron Data Library) [39] and EPDL (Evaluated Pho- 
ton Data Library) [40], as these compilations are intended to provide a consistent set of data for 
electron-photon trasport calculations. 
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